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Abstract 

Using a large sample of (j) l^^^ decays collected by the FOCUS photo- 

production experiment at Fermilab, we present new measurements of two semilep- 
tonic form factor ratios: and r2- We find r^ = 1.549 it 0.250 it 0.145 and r2 = 
0.713 lb 0.202 lb 0.266. These values are consistent with and r2 form factors mea- 
sured for the process 
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1 Introduction 



This paper provides new measurements of the parameters that describe 
(j) /i^z/ decay. The -D^ (j) fi'^u decay amplitude is described [1] by four 
form factors with an assumed (pole form) dependence. Following earlier 
experimental work [2,3,4,5,6,7,8,9,10,11,12,13], the Df — > /x+i/ amplitude 
is then described by ratios of form factors taken at = 0. The traditional 
set is: r2, r^, and which we define explicitly after Equation 1. According to 
flavor SU(3) symmetry, one expects that the the form factor ratios describing 
(f) fi^u should be close to those describing K fi^u since the 

only differs from the through the replacement of a d quark by a s 
quark spectator. The existing lattice gauge calculations [14] predict that the 
form factor ratios describing — > i'^ui should lie within 10% of those 
describing — > K i'^Ui. Although the measured form factors are quite 
consistent between (p i'^ui and K there is presently a 3.3 

a discrepancy between the r2 values measured for these two processes with 
the previously measured D'^ — >• value being a factor of about 1.8 times 
larger than the r2 value measured for K [7]. One quark model 

calculation [15] offers a possible explanation for the apparent inconsistency in 
the r2 values measured for — > K i~^u^ and Df i~^i^e- 

Five kinematic variables that uniquely describe K~ ^'^ v decay are 

illustrated in Figure 1. These are the K~ invariant mass {itlk+k-) , the 
square of the fiu mass (g^), and three decay angles: the angle between the 
and the D'^ direction in the K~K^ rest frame (6'v), the angle between the 
V and the direction in the rest frame and the acoplanarity angle 
between the two decay planes (x)- These angular conventions on 6'^ and 6'v 
apply to both the and D~ . The sense of the acoplanarity variable is defined 
via a cross product expression of the form: (P^ x P^) x {Pk- x Pk+) ■ Pk~k+ 
where all momentum vectors are in the rest frame. Since this expression 
involves five momentum vectors, as one goes from — > D~ one must change 
X — > —X ill Equation 1 to get the same intensity for the and D~ assuming 
CP symmetry. 




Fig. 1. Definition of kinematic variables. 



3 



Using the notation of [1], we write the decay distribution for — > ji^i/ in 
terms of the four hehcity basis form factors: i7+ , Hq , , Ht- 



drriKTr dq^ dcosOy dcos6e dx ' 

sin sin 9yB 



1 + cos^^) smeye'^B4,H+ 
- d - cos^^) sm9ye~'^B^H_ 



H — - 



— 2 sin 6*^ (cos 9yB^ + 

o... o... 9yB^ {e'^H+ + e-'^H_ 
+ 2 cos 9e (cos 9yB^ + Ae'^) Hq 
+ 2 (cos9yB^ + Ae'^) Ht 



2\ 



(1) 



where K is the momentum of the K~K^ system in the rest frame of the . 
The first term gives the intensity for the yU"*" to be right-handed, while the 
(highly suppressed) second term gives the intensity for it to be left-handed. 
The helicity basis form factors are given by: 



Ht{q') 



MdK 



{Ml - - q^){MD + mK^)A,{q^) - 4 



Md + rriKn 



{MD + MK^)A,{q')- 



Md + Mk^ 



-A,{q') -t- 



A2{q') 
2g2 



Md + Mk-, 



A,{q^) 



The vector and axial form factors are generally parameterized by a pole dom- 
inance form: 



Ar{q' 



Am 



1 - q^/Ml 



V{q') 



V{Q) 



1 - qVMl 



where we use nominal (spectroscopic) pole masses of Ma = 2.5 GeV/c^ and 
My — 2.1 GeV/c^. The denotes the Breit-Wigner amplitude describing the 
(f) resonance: ^ 



B. 



-ml + imoT 



^ We are using a p-wave Breit-Wigner form with a width proportional to the cube 
of the kaon momentum in the kaon-kaon rest frame (P*) over the value of this 
momentum when the kaon-kaon mass equals the resonant mass (Pq)- The squared 
modulus of our Breit-Wigner form will have an effective P*^ dependence in the 
numerator as well. Two powers P* come explicitly from the P* in the numerator of 
the amplitude and one power arises from the 4 body phase space. 
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Equation 1 includes a possible s-wave amplitude coupling to the virtual 
with the same dependence as that of the Ho (or Ht) form factor. Evidence for 
such an s-wave amplitude term for the decay —>■ K'tt^ij,'^!/ was presented 
in Reference [3] . An explicit search was made for s-wave amplitude interference 
with the process Df — > n'^v and no evidence for this interference was 
seen. We were able to limit the s-wave contribution to be less than 5% of 
the maximum of the Breit-Wigner peak in the Hq piece of Equation 1 at 
the 90% confidence level. The results presented here will therefore assume 
A — in Equation 1. Under these assumptions, the decay intensity is then 
parameterized by the = V{0)/Ai{0) , ra = ^(0)/Ai(0) , rg = ^3(0)/^i(0) 
form factor ratios describing the amplitude. Throughout this 

paper, unless explicitly stated otherwise, the charge conjugate is also implied 
when a decay mode of a specific charge is stated. 



2 Experimental and analysis details 

The data for this paper were collected in the Wideband photoproduction ex- 
periment FOCUS during the Fermilab 1996-1997 fixed-target run. In FOCUS, 
a forward multi-particle spectrometer is used to measure the interactions of 
high energy photons on a segmented BeO target. The FOCUS detector is a 
large aperture, fixed-target spectrometer with excellent vertexing and particle 
identification. Most of the FOCUS experiment and analysis techniques have 
been described previously [3,16,17,18,20]. Our analysis cuts were chosen to 
give reasonably uniform acceptance over the five kinematic decay variables, 
while maintaining a strong rejection of backgrounds. To suppress background 
from the re-interaction of particles in the target region which can mimic a de- 
cay vertex, we required that the charm secondary vertex was located at least 
three standard deviation outside of all solid material including our target and 
target microstrip system. 

To isolate the — > ii'^u topology, we required that candidate muon, pion, 
and kaon tracks appeared in a secondary vertex with a confidence level ex- 
ceeding 1%. The muon track, when extrapolated to the shielded muon arrays, 
was required to match muon hits with a confidence level exceeding 5%. The 
kaon was required to have a Cerenkov light pattern more consistent with that 
for a kaon than that for a pion by 1 unit of log hkelihood [18]. To further 
reduce non-charm background we required that our primary vertex consisted 
of at least two charged tracks. To further reduce muon misidentification, a 
muon candidate was allowed to have at most one missing hit in the 6 planes 
comprising our inner muon system and an energy exceeding 10 GeV. In order 
to suppress muons from pions and kaons decaying within our apparatus, we 
required that each muon candidate had a confidence level exceeding 1% to 
the hypothesis that it had a consistent trajectory through our two analysis 
magnets. 
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Non-charm and random combinatoric backgrounds were reduced by requir- 
ing both a detachment between the vertex containing the K~K~^fi^ and the 
primary production vertex of at least 5 standard deviations. 

Possible background from K^K^tt^, where a pion is misidentified 

as a muon, was reduced by treating the muon as a pion and requiring the 
reconstructed KKtt mass be less than 1.8 GeV/c^. The mx+K- distribution 
for our — ^ K^K~ fi^u candidates is shown in Figure 2. 

200 
150 



50 


1.00 1.05 1.10 

M(K+K") (GeV/c^) 

Fig. 2. The data is the solid histogram and ccbar background Monte Carlo is the 
dashed histogram. The ccbar background Monte Carlo is normalized to the same 
number of events in the sideband region 1.040 GeV/c^ < rrix+K- < 1.14: GeV/c^. 

The technique used to reconstruct the neutrino momentum through the 
line-of-flight, and tests of our ability to simulate the resolution on kinematic 
variables that rely on the neutrino momentum are described in Reference [3]. 




3 Fitting Technique 



The and r2 form factors were fit to the probability density function described 
by four fitted kinematic variables (g^, cos6'v, cos6'^, and x) fo^' decays in the 
mass range 1.010 < rriK+K- < 1-030. We use a variant of the continuous 
fitting technique developed by the E691 Collaboration [19] for fitting decay 
intensities where several of the kinematic variables have very poor resolution 
such as the four variables that rely on reconstructed neutrino kinematics. 

The fit which determines the r„ and r2 form factor ratios minimizes the sum 
of W7 = —2 In / where I is the normalized decay intensity at each datum. The 
intensity at each datum is estimated by the weight of Monte Carlo events 
that lie within a small window of each of the four (reconstructed) kinematic 
variables of the given datum. ^ Some care is needed in choosing a reasonable 
size of the windows since too small a window will result in fiuctuations due 

^ A Monte Carlo event had to lie within 0.08 to each datum in both cos^v and 
cos0£, within 0.18 radians in Xi arid within 0.8 GeV^/c^ in . 
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to finite Monte Carlo statistics, and too large a window will create a bias in 
the result. Our principal tool in deciding a reasonable window was to check 
for biases and fluctuations outside of reported statistical errors using a Monte 
Carlo simulation that included charm backgrounds with and r2 values very 
close to the result reported here. Variations in the final results due to the 
window choice were included in the estimate of the systematic error. 

Two Monte Carlos were employed: a weighted — > (p n^v signal Monte 
Carlo that was generated flat in the four fitted kinematic variables {q^, cos ^v, 
cos9i, and x) ^-iid an unweighted background Monte Carlo which simulated 
all known charm decays (apart from the ^ (f) fi^u signal) as well as 
our misidentification levels. The intensity about each datum is an appropriate 
average of the weighted signal MC and the unweighted background Monte 
Carlo. The averaging depended on the background fraction determined by 
matching the number of events in the (j) sideband region 1.040 < itik+k- < 
1.14 GeV/c^ in the background Monte Carlo to the sideband level observed 
in the data. 

Our fit was to the and r2 form factor ratios with and the possible s- 
wave amplitude set to zero. We decided not to fit for the form factor ratio 
since our anticipated error given our sample size would be too large to be 
meaningful. 

Figure 3 compares the data and model for projections of cos^v; cos^^, x 
^^/O'max- Most of thcsc distributions follow the predicted values reasonably 
well. A slight discrepancy is evident in the low q^/q^^^ projection (below 
q'^/q^ax < 0-2 )• A stronger effect was observed in Reference [2] possibly owing 
to our much larger yield in the K't^'^h'^v final state. 

Figure 4 compares the cos^^v and cos^^£ distribution between the data and our 
Monte Carlo model for events at high and low (f. As q^ — > g^max one expects an 
isotropic distribution in both cos and cos d(, since all three helicity basis form 
factors become equal. The data presented in Figure 4 match this expectation 
relatively well. 



4 Form Factor Ratio Systematic Errors 

Three basic approaches were used to determine the systematic error on the 
form factor ratios. In the first approach, we measured the stability of the 
branching ratio with respect to variations in analysis cuts designed to suppress 
backgrounds. In these studies we varied cuts such as the detachment criteria, 
particle identification cuts, vertex isolation cuts, and the out-of-material cut. 
Sixteen such cut sets were considered. In the second approach, we split our 
sample according to a variety of criteria deemed relevant to our acceptance, 
production, and decay models and estimated a systematic based on the con- 



7 




Fig. 3. The data are given by the upper histogram. The model (crosses with flats) 
includes the signal computed with the fitted form factor ratios plus the sideband 
normalized ccbar background. The lower histogram (dashed) shows the projections 
of the ccbar background, (a) cos^v (b) cosOi (c) q^/q^^x (d) X 

sistency of the form factor ratio measurements among the split samples. These 
included computing separate form factors for particles and antiparticles, and 
comparing the values obtained over the full data set to values obtained from 
a subset (2/3) of the data in which the target silicon [20] was operational. In 
the third approach we checked the stability of the branching fraction as we 
varied specific parameters in our Monte Carlo model and fitting procedure. 
These included varying the level of the charm background and in the size of 
the kinematic windows used to associate Monte Carlo events with data points 
in the computation of the fitting intensity. The only non-neglible systematic 
error was that due to variation in the results over the sixteen cut selections. 
Combining all three non-zero systematic error estimates in quadrature we find 
= 1.549 ± 0.250 ± 0.145 and ra = 0.713 ± 0.202 ± 0.266. 



5 Summary 



In Table 1 we compare our results to other experiments. Our weighted average 
of all the experimental results is r„ = 1.678±0.213 and r2 = 1.292±0.194 where 
systematic errors have been included. We obtain a confidence level of 44.3% 
that all 5 experiments have a consistent and a confidence level of 21.5 % 
that all r2 measurements are consistent. Figure 5 is a graphical representation 
of Table 1. 
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Fig. 4. Various cos^v and cos On projections. The data are given by the upper his- 
togram. The model (crosses with flats) includes the signal computed with the fitted 
form factor ratios plus the sideband normalized ccbar background. The lower his- 
togram (dashed) shows the projections of the ccbar background, (a) The cos 9\ dis- 
tribution for g^/g^max < 0.5 (b) The cosOy distribution for q'^/g^max > 0.5 (c) The 
cos 6*^ distribution for q^/q^m&x < 0.5 (d) The cos9e distribution for g^/g^max > 0.5 



Table 1 

Measurements of the Df (j) ^'^v form factor ratios 



Group 




r2 


This work 


1.549 ±0.250 ±0.145 


0.713 ±0.202 ±0.266 


E791[7] 


2.27 ±0.35 ±0.22 


1.570 ±0.250 ±0.190 


CLE0[8] 


0.9 ±0.6 ±0.3 


1.400 ± 0.500 ± 0.300 


E653[12] 


2.3 ± 1.0 ±0.4 


2.100 ±0.550 ±0.200 


E687[10] 


1.8 ±0.9 ±0.2 


1.100 ±0.800 ±0.100 



Our measured and r2 values for n^u are very consistent with 

our measured and r2 values for ^ K*^ jX^v [2]. The measurements 
reported here call into question the apparent inconsistency between r2 values 
the K i'^i'i and form factors present in previously 

published data and are consistent with the theoretical expectation that the 
form factors for the two processes should be very similar. 
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Fig. 5. Comparision to previous data and calculations of the Df — > (p fi^v form 
factors. The three solid green lines represent the weighted average of the world's 
experimental data for (p and ±1 o". The three dashed blue lines represent 

the weighted average of the world's experimental data on K ^ru and ±1 a. 
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